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Abstract 

We revisit inflationary cosmology of axion models in the light of recent developments on the 
inflaton decay in supergravity. We find that all the cosmological difficulties, including gravitino, 
axino overproduction and axionic isocurvature fluctuation, can be avoided if the saxion field has 
large initial amplitude during inflation and decays before big-bang nucleosynthesis. 
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I. INTRODUCTION 



Although the standard model of the particle physics has achieved great successes, it 
still has some theoretical problems. One is the gauge hierarchy problem, and another is 
the strong CP problem. On the otherhand, cosmological problems such as dark matter and 
baryon asymmetry of the universe cannot be explained within the framework of the standard 
model. These problems indicate that there is an underlying physics beyond the standard 
model. 

Supersymmetry (SUSY) is one of the best motivated candidates as physics beyond the 
standard model to solve the gauge hierarchy problem. On the other hand, Peccei-Quinn (PQ) 
mechanism is a simple promising solution to the strong CP problem As a consequence 
of the breaking of PQ symmetry, the existence of a pseudo-Nambu-Goldstone boson, axion, 
is predicted [2J. Thus we are eager to study the axion models within the framework of 
SUSY. However, cosmology of SUSY axion models is highly non-trivial. It is well known 
that gravitino, which is the fermionic superpartner of the graviton and has a long lifetime, 
is produced during reheating after inflation and its decay gives serious effect to thermal 
histry of the universe (this is called gravitino problem). In addition to the gravitino, there 
appear other long-lived particles called saxion, the scalar partner of the axion, and axino, the 
fermionic superpartner of the axion [3j. Both of them have a potential to cause cosmological 
disaster, which may be even more problematic than the usual gravitino problem 



Furthermore, recently it ispointed out that gravitinos are also produced non-thermally 
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from the inflaton decay [10|, llll, ll2|, ll3|, LU, l!5| . Inclusion of such contributions makes the 
gravitino problem much worse so that some inflation models may be excluded depending on 
the mass of the gravitino. 

This is not the end of the story. The axion field induces an isocurvature fluctuation 16] , 
whose amplitude is proportional to the Hubble scale during inflation. Recent cosmological 
observations are consistent with pure adiabatic one as the initial density fluctuation of 



the universe 



171 ]. and hence the fraction of the isocurvature component to the density 



perturbation is constrained. This excludes the high-scale inflation models such as chaotic 
inflation or many hybrid inflation models. 

Thus SUSY axion models seem to suffer from various cosmological problems, which ex- 
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elude many inflation models. However, in this paper we point out that a resolution to all the 
above difficulties is already built in SUSY axion models themselves. The saxion is naturally 
expected to have a large initial amplitude of order of the reduced Planck scale M P during 
inflation, and starts coherent oscillation after inflation. Soon after reheating due to the 
inflaton decay, the saxion dominates the universe since the energy density of the saxion is 
almost comparable to the total energy density during the inflaton-dominated era. Finally, 
the saxion decays at later epoch and reheats the universe again releasing huge entropy. A 
remarkable feature is that huge entropy produced by the saxion decay dilutes the possibly 
harmful gravitinos and axinos. This provides a solution to the cosmological gravitino and 
axino problems. Moreover, the amplitude of the axionic isocurvature fluctuation is also sig- 
nificantly reduced for a large initial amplitude of the saxion, which makes high-scale inflation 
models such as chaotic inflation compatible with cosmological constraints. 

This paper is organized as follows. In SecJTT]we overview cosmological problems with the 
gravitino, axino and axion. In Sec lHIl we show that entropy-production by the saxion decay 
overcomes these cosmological difficulties. We conclude in SecJIVi 



II. COSMOLOGICAL PROBLEMS 



In this section we briefly summarize cosmological constraints on the abundances of long- 
lived particles, the gravitino, axino and axion, which appear in SUSY axion models. 



A. Gravitino 



Gravitinos are produced in the early universe through scatterings of particles in thermal 



bath. For an unstable gravitino, its decay may significantly affect Big-Bang nucleosynthesis 
(BBN) [19I . [20! . If the gravitino is stable, it may have too large contribution to the present 
matter density of the universe [2l|]. Thus the gravitino abundance is constrained for nearly 



all mass range. The abundance is calculated as 20 
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where rrig is the gluino mass and Tr is the reheating temperature of the universe defined as 
T R = (10/7r 2 ^) 1 /4 V r totajM p. 



Recently, it was pointed out that gravitinos are a 



the inflaton through supergravity effects 



3, 
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so produced directly from the decay of 
, LL3[ . Taking account of the mixing of 



the inflaton with SUSY breaking field, the decay rate of the inflaton into gravitino pair is 
written as 

where ttu denotes the inflaton mass and Q§ is the effective coupling to the gravitino given 
in Ref. 12j. In this paper we assume a dynamical SUSY breaking scenario with dynamical 
scale A. If <C mz where mz denotes the mass of the SUSY breaking field (hereafter for 
simplicity we assume m z ~ A), the direct production process cannot be suppressed. In this 
case, the effective coupling is given by \Q^\ ~ 3((0)/Mp)(m 3 / 2 /"T.^) and the decay rate is 
estimated as [3, Q 

This gives the gravitino abundance as 
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If > A, the direct decay of the inflaton into gravitino pair can be suppressed if there 

are no couplings such as 5K ~ \<fi\ 2 zz in the Kahler potential. Instead, the inflaton decays 

into hidden gauge sector through the anomaly effects, and each hidden hadron eventually 

produces at least one gravitino [14j . The partial decay rate of the inflaton into hidden gauge 

sector is given by 

r , ~^(T h -T h ?(^-\ 2 ^ (5) 

anomaly - 35^3 * R ) ^ j M y W 

where N h is the number of generators and is the gauge coupling constant of the hidden 

gauge group. Tq and are the Dynkin index of the adjoint representation and matter 

fields in the representation of the dimension d R , respectively. Then the abundance of non- 

thermally produced gravitinos is given by 

3Tp 



where e is 0(1) constant given by e = N 3/2 N h a 2 h (T£ - T£) 2 (here N 3/2 denotes the averaged 
number of produced gravitinos per hidden hadron jet). 

If the gravitino is unstable, photons or hadrons produced through the decay of the grav- 
itino may affect light element abundances synthesized through BBN. On the other hand, if 
the gravitino is stable, it contributes to the present dark matter density. Both set the upper 
bound on the gravitino abundance and reheating temperature of the universe. Moreover, we 
can see that is proportional to Tr, while iv^ * s proportional to T^ 1 . Thus lowering 

Tr does not ameliorate the situation. Also it was found that the reheating temperature of 
the universe is bounded from below, due to the spontaneous decay processes of the inflaton 



through top Yukawa coupling 23], as 



/ %o TV , f 228.75 m v yfi / (0) 

where yt is the top Yukawa coupling. Including non-thermally produced contribution, many 
inflation models are severely constrained. For example, if the gravitino mass is ~ 1 TeV 
and its hadronic branching ratio is 0(1), typical inflation models such as new, hybrid and 



chaotic inflation models are excluded 15] . 



B. Axino 



In SUSY axion model, the axion forms a supemultiplet which includes scalar partner of 



the axion called saxion, and fermionic superpartner of the axion cal 



ed axino. Both of them 
Similar to the 
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have long lifetime and may significantly affect cosmology 
gravitino, axinos are also produced through scatterings of the particles in thermal bath. The 
resulting abundance is calculated as \^ 

where g s is the QCD gauge coupling constant. This expression is valid for Tr > 10 TeV, 
and the abundance is suppressed for T R < 1 TeV because SUSY particles are not produced 



efficiently. Althoug 



i the axino mass is model dependent, it may have the mass of order of 



the gravitino mass 24j . Here we assume that the axino mass is of the order of the gravitino 
mass. 
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From (jSj), one can see that very low reheating temperature is needed in order to avoid the 
overproduction of the axino if it is the LSP. This sets an upper bound on Tr as Tr < 1 TeV 
for m- a > 1 GeV, and T R < lTeV(lGeV/m s ) for m~ a < 1 GeV when F a = 10 11 GeV. 



Axion 



The axion is a pseudo-Nambu-Goldstone boson associated with spontaneous breaking 
of the PQ symmetry. The axion is practically massless for T > 1 GeV due to the finite- 
temperature effect, and begins to oscillate after it becomes massive for T < 1 GeV. The 
abundance of axion in the form of coherent oscillation is estimated as 



25 1 
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where 9% denotes the initial misalignment angle of the axion, which is naturally expected 
to be 0(1) without fine-tuning. Thus the upper bound on the PQ scale is F a < 10 12 GeV, 



although late-time entropy production can relax this upper bound 26| . On the other hand, 



astrophysical arguments require that F a should be larger than ~ 10 10 GeV (2 

If the PQ symmetry is broken during or before the inflation, the quantum fluctuation of 
the axion induces isocurvature fluctuation with magnitude ~ Hi/ {itSi) [if]], where s« denotes 
the field value of the saxion during inflation. Recent cosmological observations indicate that 



the ratio o 
about 0.3 



' the magnitude of isocurvature perturbation to adiabatic one should be less than 
17,0]. This leads to the constraint on the Hubble scale during inflation, 



0.13 J \F a J \10 12 GeV^ 

Thus, high-scale inflation models are excluded unless Sj is as large as the Planck scale. As 
a result, low-scale inflation models such as a new inflation model were considered to be 
suitable for the axion cosmology. Note that if Hj > F a , the PQ symmetry may be restored 
during inflation and hence isocurvature constraint does not apply, although there may be 
problematic domain wall formation for generic axion models. 

In Fig. [T]the axion isocurvature constraint is shown on the m^-(4>) plane with Sj = F a = 
10 12 GeV. 1 Also we show the prediction of typical inflation models in supergravity, new 



Here we have estimated the inflation energy scale by Hi ~ m^(tp) /y/3Mp. For inflation models with the 
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FIG. 1: Axion isocurvature constraint on the mj,-{(j>) plane with Sj = F a = 10 12 GeV. The shaded 
region is excluded from the isocurvature perturbation constraint. Green solid lines represent new 
inflation models with n = 4, 5, 6 from left to right, red dotted line represents chaotic inflation model, 
purple dot-dashed line represents hybrid-inflation model and blue dashed lines represent smooth- 
hybrid inflation models for m = 2-6 from left to right (for definitions of n and m see footnote). 
Region above thin solid black lines are excuded from gravitino overproduction in inflaton decay for 
(a) m 3/2 = 1 TeV and (b) m 3/2 = 100 GeV, with T R = 1 TeV. 



28|, [29], hybrid 30|, l3l|] , smooth-hybrid 32], and chaotic inflation [33] models. 2 Furthermore, 
non-thermally produced gravitinos restrict the parameter region. Regions above thin solid 
black lines are excluded from gravitino overproduction, for (a) m 3 /2 = 1 TeV and (b) 
m 3/2 = 100 GeV. Here we set Tr — 1 TeV, since axino thermal production sets the upper 
bound on the reheating temperature as Tr < 1 TeV. Saxions are not harmful for cosmology 
for these parameter sets. We can see that most inflation models are excluded. However, if one 
allows a large field value of the saxion during inflation, the effective PQ scale during inflation 



potential like V ~ (v 2 — g4>") 2 , as is the case for new, hybrid and smooth-hybrid inflation models, this 
estimation is correct except for the numerical factor n. For the chaotic inflation model, this evaluation 
cannot be applied and the correct estimate is Hj ~ 10 14 GeV. 
2 We assume the following superpotentials for new and smooth hybrid inflation models : W = %jj(v 2 — gcf> n ) 
for new inflation and W = ip[v 2 — (4><j>) m / M 2,n ~ 2 ] for smooth hybrid new inflation model. 
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can be practically much larger than F a , which suppress the magnitude of the isocurvature 
perturbation 34j, as we will see in the next section. 



III. SAVING INFLATION MODELS 

All the above cosmological difficulties are avoided by taking account of the dynamics of 
the saxion. In the above arguments, we have assumed that there was no entropy production 
after reheating by the inflaton. However, in SUSY axion model, the saxion can have the large 
initial amplitude of order of the Planck scale. For example, let us consider the case where 
the saxion has negative Hubble mass term. Since the saxion s corresponds to a flat direction 
in the scalar potential, the saxion field rolls away due to the negative Hubble mass term 
during inflation until the field value becomes the Planck scale where the potential becomes 
steep as ~ exp(s 2 /Mj,). Such saxion condensate provides a source of late-time entropy 
production [35|] and dilutes the harmful gravitino and axino. Furthermore, large initial 
amplitude during inflation suppresses the axionic isocurvature fluctuation. Therefore, there 
arises a possibility that many inflation models are free from cosmological difficulties which 
we encountered in the previous section. 

A. Saxion with large initial amplitude 

The saxion corresponds to a flat direction along which the PQ scalars do not feel the 
potential. A flat direction is lifted by the SUSY breaking effect and hence the saxion mass 
(m s ) is naturally expected to be of order of the gravitino mass. 

The coherent saxion oscillation starts at H ~ m s , with the initial amplitude If 
Si ~ F a and the saxion is relatively light, its decay causes various cosmological difficulties, 
which severely restricts the saxion abundance (si, El- As a result, the upper bond on the 
reheating temperature becomes much more stringent than the case where the saxion and 
axino are absent for wide range of the saxion mass [9|. 

Here we consider another possibility. Thermal history of the universe can be significantly 
modified due to the saxion coherent oscillation with a large initial amplitude ~ Mp. In 
this case the saxion dominates the universe immediately after the inflaton decays, and finally 
the universe is reheated again by the saxion decay. However, if the saxion main decay mode 



S 



is into two axions (s — > 2a), the produced axions contribute to the extra relativistic degrees 
of freedom, which changes the Hubble expansion rate and the BBN prediction. Thus we 
need to investigate the saxion decay modes in order to ensure that the saxion decay does 



not produce too many axions. The decay rate into two axions is estimated as 



24] 



f ^ 7773 

where / = ^2 { qfvf / F% with the VEV of the z-th PQ scalar field V{ and its PQ charge 
This often gives the dominant contribution to the total saxion decay rate if / ~ 1. In order 
to realize consistent cosmology, this must not be the main decay mode. 

Another important decay mode is into gluons, arising from the coupling through the QCD 
anomaly effect. The decay rate is estimated as 

at mi 



r.(.- w ) = 3^- (12) 



For the DFSZ axion model 



361 ]. the saxion can also decay into fermion-anti-fermion pair 
through the coupling arising from mixing of the PQ scalar with MSSM Higgs doublets. The 
decay rate is estimated as 

o \2 / \ 2 / a 2 \ 3/2 

2x \ ( m di \ ( , 47774 
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for the decay into down-type quarks di (i = 1, 2, 3) where x = tan (3 = (H u ) / (Ha), and 



7r \x + x 1 / V F a J \ m 



for the decay into up-type quarks (i = 1,2,3). The decay rate into Higgs boson pair is 
also comparable, 

1 3 / \ 4 / 4 2 \ 1/2 

1 / /i \ / 4m^ x 



where /x-parameter is of order of the weak scale and set to 300 GeV here. For simplicity we 
consider only the decay into the lightest Higgs bosons. Decays into other Higgs bosons also 
have comparable rate if kinematically allowed. Also we neglect decay into SUSY particle 
pair assuming that such decay modes are kinematically forbidden. 

In order not to contradict with observations, the increase of the effective number of 
neutrino species AN U should be smaller than about 1, which restricts the branching ratio 



into axions B a as 
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FIG. 2: Branching ratio into axions in DFSZ model for / = 1 (upper) and / = 0.1 (lower). Region 
above the thin dashed line is excluded. 

where AiV^bound) ~ 1, and T s denotes the decay temperature of the saxion. 

In Fig. [21 B a as a function of the saxion mass with / = 1 and / = 0.1 is shown. As 
a representative value, we show B a = 0.2 by the thin dashed line above which too many 
axions are produced by the saxion decay and the constraint (fl6|) is not satisfied. It can be 
seen that for some mass ranges, the decay into axions is subdominant process. For that case 
the saxion is an ideal candidate of the source of late-time entropy production. 

B. Dilution of gravitinos and axinos 

We have shown that late-time entropy production by the saxion decay can can take 
place without producing too many axions. Now let us investigate the gravitino and axino 
abundances after the dilution by the saxion decay. 

If we assume the main decay mode is s — > bb, the typical decay temperature of the saxion 
is estimated as 



harmful axinos and gravitinos, produced either thermally or non-thermally, are diluted by 
the saxion decay. After the dilution, thermally produced gravitinos during reheating pro- 
cesses after inflation have only negligible abundance of order YSfJ < 10 -20 independently of 
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Tr. On the other hand, the gravitino abundance produced from inflaton decay is estimated 

as 
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for > A where 7 is defined as 



7 = { 1 f ° r T « =* Tfl , (20) 

Tr/T osc for T osc < Tr 

with the temperature at which the saxion starts to oscillate T osc . Axinos are also diluted by 
the saxion decay and the resultant abundance is estimated as 3 

-1.3 (, / Fa \ 2 ( T s \ ( Mp N 



In Figs. [3] and HI the allowed and excluded region in m^-l^tfi) plane are shown for the unstable 
gravitino with 777,3/2 = 1 TeV and stable gravitino with m 3 / 2 = 100 GeV. We put a constraint 
^3/2 < 10~ 16 for the former case Note that axinos decay well before BBN for rria ^ 1 TeV 
and do not affect BBN. Instead, the LSP abundance emitted from the axino puts a constraint 
on the axino abundance, and this is also satisfied. It can be seen that all of the inflation 
models are allowed for the reheating temperature Tp > 10 9 (10 6 ) GeV for m 3 / 2 = 1 TeV 
(100 GeV). Moreover, the constraint on the inflation scale ffTU]) is relaxed for s, ~ Mp, 
which enables even the chaotic inflation model to be consistent with observations. 4 



3 Here we have assumed that the reheating due to the inflaton decay completes after the saxion oscillation. 
Otherwise, the axino abundance is suppressed by sf ~ Mp, not F%, and is safely neglected. 

4 Again it should be noted that isocurvature constraint on the chaotic inflation model cannot be read from 
these figures. See the footnote in Sec lII CI 
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FIG. 3: Constraints on inflation models for m 3/ / 2 = 1 TeV and F a = 10 12 GeV with Si = Mp. 
Region above the thin solid black lines is excluded for (a) Tr = 10 9 GeV, (b) Tr = 10 7 GeV and 
(c) Tr = 10 5 GeV. The shaded region is excluded from axionic isocurvature constraint. 
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FIG. 4: Same as Fig. H except for m 3/2 = 100 GeV and F a = 10 11 GeV 
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Thus gravitinos are sufficiently diluted so that they do not cause cosmological difficulties. 
Interestingly, for the gravitino mass m 3 / 2 ~ 100 GeV, non-thermally produced gravitinos 
can have the desired abundance as dark matter of the universe. One issue to be addressed 
is the free-streaming length of the gravitino dark matter. It is calculated as 
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where R eq is the comoving Hubble scale (~ 108Mpc) and u eq = v/y/l — v 2 with the veloc- 
ity of the gravitino v, both are evaluated at the time of matter-radiation equality. u eq is 
estimated as 

1/3 / m \ -4/3 
j_ 



u cq ~ 8 x 10 
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Hence the free-streaming length is negligibly small and the gravitino dark matter acts 
as cold dark matter. The axion is diluted by the saxion decay because T s < 1 GeV for 
m s ^100 GeV, but it may still have the comparable abundance to the present dark matter 
abundance. 

For m s > 1 TeV, there arises a possibility that the saxion decays into SUSY particles 
if the decay is kinematically allowed. According to Ref. {38], the decay rate of the saxion 
into gauginos are roughly the same as that into gauge bosons. Thus in general LSPs are 
overproduced by the saxion decay. However, for the LSP with rather large annihilation cross 
section, these non-thermally produced LSP abundance is reduced. The abundance is given 
by 

2m L sp Ps 



Plsp 



mm 
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(24) 



/;.'., .s \/ 8ir 2 g^(T s ) (av)T s Mp 
where B s denotes the branching ratio of the saxion into SUSY particles, and (av) denotes 
the thermally averaged annihilation cross section of the LSP. For the LSP with (av) ~ 
10~ 7 GeV -2 , as is realized for the case of wino- or higgsino-like LSP with the mass of 
O(100) GeV 39j, or bino-like LSP in the S"-channel resonance region [40(, the resultant 



LSP abundance can account for the present dark matter of the universe. Also there is 
contribution from the axion coherent oscillation to the dark matter abundance. Thus the 
dark matter may consist of a mixture of the axions and non-thermally produced LSPs. 
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C. Baryon asymmetry 



Note that baryon asymmetry is also diluted by the saxion decay. Here we show that 



Affleck-Dine mechanism 



411 ] can create appropriate amount of baryon asymmetry. 5 Let us 
denote the Affleck-Dine (AD) field, which is one of the flat directions in the scalar potential in 
the MSSM, as if). The flat direction is lifted by SUSY breaking effects and non-renormalizable 
terms in the superpotential IUnr = ip n /nM n ~ 3 where n(> 4) is an integer and M denotes 
the cutoff scale. The resultant baryon asymmetry created by the coherent motion of the AD 
field is estimated as 

n B _ nB_ps_ ^ <?cpra3/2|^oB| 2 3r 8 
s p s s ~ rn^sj 4 

where #cp denotes the effective CP angle, which is naturally expected to be 0(1), and ^ os 
is the field value at the onset of the oscillation of the AD field. If there exists a negative 
Hubble mass term for the AD field, the field value is given by |VU ~ (m^M n - 3 ) 1/(n ~ 2) - For 
the specific case n = 6, we obtain 

n B o m-njc f^Y 3/2 f m m\f T s \fM\ 3/2 fM P \ 2 

and hence AD mechanism works well. Although relatively large Q-balls are formed through 
the AD mechanism (Q ~ 10 20 ) 2], their cosmological effects are safely neglected, since 
Q-balls decay before the saxion decays and LSPs emitted non-thermally by the Q-ball decay 
is diluted by the entropy-production from the saxion. 6 



IV. CONCLUSIONS 

We have shown that in SUSY axion models the saxion decay naturally dilutes the axino 
and gravitino abundances produced both thermally or non-thermally if the initial amplitude 
is of the order of Mp. The large initial amplitude of the saxion also saves high-scale inflation 
models such as chaotic inflation and hybrid inflation models from producing too much axionic 



5 See e.g., Ref. [42J for more details in the case of late-time entropy production. 

6 The Affleck-Dine mechanism described above has the baryonic isocurvature problem if the Hubble induced 
A-term is absent as expected for most inflation models [44(. In this case high scale infaltion models are 
disfavored. However, the problem is avoided for AD mechanism without superpotential (4^ . 
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isocurvature fluctuation. Note that the preexisting baryon asymmetry is also diluted by the 
saxion. But the Affleck-Dine mechanism can create large baryon asymmetry which survives 
the dilution. 

Finally we mention another possibility to dilute the gravitino and axino abundances. 



One may consider that thermal inflation 45j, |46| can provide sufficient dilution. If the PQ 



scalar is trapped at the origin during inflation and remains there until later epoch due to 
the finite-temperature effect, the PQ scalar itself can cause a thermal inflation for some 
class of models 47]. Axionic isocurvature fluctuation does not arise since the PQ symmetry 



is restored during inflation. However, we must take care of topological defects formation 
after the onset of the PQ scalar oscillation, which may spoil the subsequent cosmological 
evolution scenario. 
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